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ABSTRACT 

We present 19 nearby (z < 0.3) SCR radio galaxies imaged at low- and high-excitation as part of a Cycle 
15 Hubble Space Telescope snapshot survey with the Advanced Camera for Surveys. These images consist of 
exposures of the Ha (6563 A, plus [Nil] contamination) and [Olll]A5007 emission lines using narrow-band 
linear ramp filters adjusted according to the redshift of the target. To facilitate continuum subtraction, a single- 
pointing 60 s line-free exposure was taken with a medium-band filter appropriate for the target's redshift. We 
discuss the steps taken to reduce these images independently of the automated recalibration pipeline so as to use 
more recent ACS flat-field data as well as to better reject cosmic rays. We describe the method used to produce 
continuum-free (pure line-emission) images, and present these images along with qualitative descriptions of 
the narrow-line region morphologies we observe. We present Ha+[Nll] and [Olll] line fluxes from aperture 
photometry, finding the values to fall expectedly on the redshift-luminosity trend from a past HST/WFFC2 
emission line study of a larger, generally higher redshift subset of the SCR. We also find expected trends 
between emission line luminosity and total radio power, as well as a positive correlation between the size of 
the emission line region and redshift. We discuss the associated interpretation of these results, and conclude 
with a summary of future work enabled by this dataset. 

Subject headings: galaxies: active — galaxies: emission lines — radio continuum: galaxies 



1. INTRODUCTION 

Characteristically intense nuclear emission in radio galax- 
ies, along with the highly collimated jets powered by ac- 
cretion onto their central engines, can influence star forma- 
tion, excitation and ionization of the ISM, as well as pro- 
vide kinetic stresses on the hot X-ra y emitting coronal gas 
that pervades clusters of galaxies d OuiHen & Bower 1999; 
i McNamara et al.l l2000l; 'Bl anton et al.. 2001: Reynolds et al 
I2OO2; Ruszkowski & Begel manll2002|) . Narrow-band imag- 
ing surveys of radio galaxies at redshifts z > 0.6 were among 
the first studies to establish such a connection, having shown 
the optical line-emitting gas (T ^ 10^ K) near the nuclei of 
these objects to spatially alig n with the radio jet axes on kilo- 
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propagation of t he radio jet, as well as photoionization from 
the AGN (e^ lBaum & HeckmanI [T989al: iMcCarthvl [T993I: 
iBest et all l2000h . While these shocks can trigger star for- 
mation along the regions of the ISM excited into line emis- 
sion (Chambers e t al. 1987) , more recent work has suggested 
that feedback from the AGN may also play a role in ulti- 
mately quenching the star formation rate (SFR) in the late 
stages of the host galaxy's evolution, thereby ushering its 
rapid passage from the 'blue cloud' to the 'red sequence' (e.g., 
ICowie et al. 1996 : BelTet al.ll2004l: iFaber et a l. 2007. and ref- 
erences therein). Work is ongoing in reconciling the seem- 
ingly competing roles of AGN feedback in both inducing, 
then possibly truncating star formation at successive stages 
of galactic evo l ution ( e.g., |Si lk & Rees 1998"; 'Fabian' '1999"; 
I di Matte o et"aD 120051; iHopkins et al. 2005; Silverman et alj 
l2008h . Regardless, it is clear that the AGN a nd its host galaxy 
share codependent evolutionary paths (e.g. , 'Magorr ian et al.l 
119981: iFerrarese & MerrittI l2000l: iGebhardt et al' 20(©, and 
that studies of the ionized nuclear gas in radio galaxies, as 
unique probes of the early universe, will be key to understand- 
ing this relationship on size and time scales small and large. 

The Hubble Space Telescope (HST) has undertaken system- 
atic surveys of th e SCR catal og of radio galaxies in the UV 
('Allen etal."2002), optical (MarteLeLalJ|1999|; Ide Koff et al " 
| 2000 : P rivon et al. 2008), and near-infra red ^Madrid et al. 
p006l: iDonzelli et al.lf2007l: [Tremblav et al.ll2007TfFlovd et al 
'2008*). These programs, in concert with a robust array of 
ground-based observations of the SCR in nearly all wave- 
length regimes, have established a uniform database of cross- 
spectrum imaging and spectroscopy for a sample that is nearly 
complete with redshift and flux limit, unbiased with respect 
to optical/IR wavelengths, and containing galaxies exhibiting 
a variety of intrinsic characteristics. The completeness and 
diversity of the sample has enabled studies into radio-loud 
unification models and the dichotomy between relatively low 
power, edge-darkened Fanaroff and Riley class I radio galax- 
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ies (hereafter FR L lFanaroff & Rileyll 19741) and higher power , 
edge-brightened FR II radio galaxies lchiaberge et al.ll20QQ[ 
l2002h . if^T observations of the 3 CR have also led t o discov- 
eries of new optical and IR jets (iFloyd et al.ll2006h . a trend 
between nuclear dusty disk inclination and host galaxy ispho- 
tal shapes (Tre mblav et al. 2007) and face on disks with opti- 
cal jets ( Spark s et al.ll2000l) . Jet/disk orien tations were studied 
at unprecedented sp atial resolution (e.g., Schmitt et al.'l 2002l: 
IVerdoes Kleijn & d e Zeeuw 2005; Tremblay et al. 2006i and 
references therein). Allen et al. (2002) also found complex 
UV morphologies unlike those seen in quiescent hosts, in- 
dicative of ongoing star formation. 

In this paper we present narrow-band emission line obser- 
vations of 19 nearby (z < 0.3) SCR radio galaxies imaged at 
low- and high- excitation with the Advanced Camera for Sur- 
veys (ACS) aboard HST. These images of the warm optical 
line (Ha and [Olll]A5007) emitting gas present in the nuclei 
of the radio galaxies in our sample offer the highest sensitivity 
and spatial resolution yet available in the HST SCR database. 
This enables detailed studies of shocked and star-forming re- 
gions of the ISM in relation to radio jets, providing the already 
rich SCR dataset with a larger framework from which to study 
the various phenomena discussed above. 

We organize this paper as follows. In §2 we describe these 
HST observations in detail, and in §S we discuss the asso- 
ciated data reduction. In §4 we present the Ha+[Nll] (low- 
excitation) and [Olll] (high-excitation) emission line images, 
as well as qualitative descriptions of each. We also provide 
emission line luminosity results and an associated discussion. 
We summarize this work and discuss future studies enabled by 
this new dataset in §5. Throughout this paper we use Hq = 1\ 
km s-i Mpc-\ VLm = 0.27, and = 0.7S. 

2. SAMPLE SELECTION & OBSERVATIONS 

In this paper we present 19 radio galaxies from the low- 
redshift (z < 03) extra^alac tic subset of th e Revised Third 
Cambridge Catalog (SCR, iBennettI [l962blial : ISpinrad et al.1 
fT985h . SCR sources are selected by radio flux density at 178 
MHz, at which extended, unbeamed lobes are detected irre- 
spective of the radio jet's orientation with respect to the line of 
sight. The sample is therefore free from orientation bias and is 
nearly complete with redshift, providing a strong framework 
for statistical analysis important for studies of (e.g.) radio- 
loud unification models. Moreover, the sample consists of 
radio galaxies that exhibit a wide variety of characteristics 
(size, shape, jet morphology, radio luminosity, FR type, dust 
content, etc.), allowing for physical comparisons to be made 
among properties intrinsic to various species of radio-loud 
AGN. We select the low-redshift subset of the SCR so as to 
achieve the highest possible spatial resolution, in addition to 
requiring the Ha line to be within the redshift range of high 
sensitivity to provide uniform line ratio mapping. 

Of the 1 16 extragalactic SCR sources with z < O.S, 98 were 
awarded for observation in the Cycle 15 ACS snapshot pro- 
gram 10882 (PI: Sparks). The remaining 18 targets that were 
not requested in our proposal were excluded for various rea- 
sons, mainly due to emission line data already being present 
in the archive for that object. Ultimately, 20 of our awarded 
targets were observed prior to failure of ACS side two elec- 
tronics in January 2007. Of these observations, 19 were suc- 
cessful. Our observation of SC S71 was carried out but un- 
successful due to failure of guide star acquisition, resulting 
in unusable data. The sample consists of 4 FR I radio galax- 
ies, 12 FR II radio galaxies, and S targets exhibiting compact 
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□ Underlying z < 0.3 SCR sample. 
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Fig. 1. — Redshift distribution of the low redshift (z < 0.3) extragalactic 
subset of the SCR sample (116 galaxies), compared with the distribution of 
targets successfully imaged (prior to ACS failure) as part of the ACS snap- 
shot program presented in this paper (solid, 19 galaxies). While this dataset 
represents a small fraction of the underlying sample, it spans nearly the entire 
range of redshift. Collectively, the vast majority of the underlying sample has 
been imaged as part of past HST programs in the UV, optical, and near IR. 
See §1 for a summary of these //^r programs. 

Steep spectrum (CSS) or undiscernible radio morphology. 

While only a small subset of our initially proposed sample 
was observed, the nature of HSTs snapshot mode means that 
observations were scheduled by convenience with respect to 
HSTs other programs at the time, and not by any selection 
effect. Our group of galaxies is therefore as unbiased as the 
SCR itself. Moreover, the 19 galaxies successfully observed 
span nearly the whole redshift range (0.0075 < z< 0.224, see 
Fig. [T]) of the underlying sample and, as we will discuss in 
more detail in §4, includes galaxies exhibiting a wide array of 
emission line characteristics. 

Our 19 successful observations consist of ACS Wide Field 
Channel (WFC) exposures of the Ha+[Nii] and [Oiii]A5007 
emission lines using narrow-band (2% bandpass) linear ramp 
filters (LRFs) at wavelengths adjusted according to the red- 
shift of the target (see Table [B. The LRFs aboard ACS fea- 
ture effective central bandwidths that vary as the filter is ro- 
tated across the FOV. While the LRFs are useful in providing 
the observer narrow-band imaging capabilities at a wide va- 
riety of central wavelengths, there is a drawback in that the 
target is imaged monochromatically only over a small region 
(^ 40'' X 80'0 of the instrument's FOV. In planning the ob- 
servations for this dataset, we ensured that the target galaxy 
nuclei always corresponded to the very center of this region. 

With each target fully observed in one orbit, an exposure 
time of 2 X 200 s and 2 x 250 s was used for the Ha and [Oiii] 
lines, respectively, with a two point dither pattern (of line 
spacing 0'M45) to aide cosmic ray and hot pixel rejection. 
So as to facilitate continuum subtraction from the emission 
line images, a 60 s, single pointing exposure was taken with 
a medium (9%) bandpass ramp filter centered at rest-frame 
5500 A, covering a region of the optical spectrum absent of 
significant contamination from line emission. Subarrays were 
used to increase the efficiency of these observations, allowing 
4 images to be stored in the ACS buffer at once. The WFC 
on ACS consists of two 2048 x 4096 pixel CCDs (WFCl and 
WFC2), each with a plate scale of 0''.049 per pixel and a field 
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Fig. 2. — A representative example of our filter selection strategy. We have convolved template LINER (chosen only for example Ho; and [Nil] lines) and QSO 
(for [Oiii]) spectra through the filters chosen for our observations of 3C 234 (z = 0.18) as an example. The spectra are used only as an illustration, and are not 
meant for quantitative discussion. The counts scale is therefore arbitrary. See §2 for details on our filter selection strategy. 



of view (FOV) of 202'' x 202'' (further details may be found 
in the ACS Instrument Handbook^ See Table [T] for an ob- 
servation log related to this program. 

It is important to note that, for all "Ha imaging" presented 
in this paper, the bandpass used was contaminated by [Nil] 
emission at 6548 and 6583 A. References to Ha will there- 
fore read "Ha+[Nll] " throughout this paper. Though it is 
beyond the scope of this work to subtract this contamination, 
[Nil] flux does not impact the results presented in this paper 
in any way as we present all corellations with respect to 
Ha+[Nll] and never Ha alone. In follow-up papers providing 
analysis of these data, howe ver, new ground-based o ptical 
spectroscopy of the SCR by Buttiglione et al.' ('2009') may 
be used to better quantify the contribution of [Nil] through 
each Ha bandpass. The Buttiglione spectra, obtained using 
the Telescopio Nazionale Galileo (TNG), may be convolved 
with the spectral response of the HST ACS ramp filters, 
providing a constraint on the relative strength of Ha and 
[Nil] in the filter throughput. In doing this, one must take 
into account the fact that the TNG spectra are extracted from 
a 2"x 2"region, while many of our galaxies (e.g., 3C 33) 
are far more exten ded. In Fig. [3] we pl ot the de-reddened 
flux ratios from Buttig lione et al. I (l2009l) of the diagnostic 
line [Nll]A6584 to Ha, vs. redshift, for each galaxy in our 
sample (with the excep tion of 3C 132). The results from 
iButtiglione et al.l (12009 ). some of which are apparent in 
Fig. [51 show that the relative contributions of [Nil] and Ha 
vary widely from source-to-source, and that [Nil] appears to 
dominate in most high-excitation galaxies (HEG, see Table[2]). 



11 Boffi, F.R., et al. 2007, "ACS Instrument Handbook", Version 8.0, (Bal- 
timore: STScI). 



3. DATA REDUCTION 

Here we describe the steps taken to produce data that has 
been calibrated, cleaned of cosmic rays (CRs) and hot pixels, 
and finally subtracted of contamination from continuum emis- 
sion. The if On-the-Fly Recalibration (OTFR) pipeline uti- 
lizes pre-launch (ground) flats in the reduction of ACS LRF 
images, the use of which can contribute upwards of 7% uncer- 
tainty to reduced data given the evolution of the instrument in 
the years since launch. By reducing each image "manually", 
we are able to avoid introducing this uncertainty, as well as 
implement a more proactive approach to cleaning the images 
of CRs and hot pixels, particularly in the case of the single- 
pointing continuum images for which automated CR rejection 
via combination of dithered images is not possible. For these 
reasons the images presented in this paper have been reduced 
independently of the OTFR pipeline, in a process we describe 
below. 



3.1. Calibration 

The calibration stage begins with the raw exposures and 
corrects for flat-fielding, bias and dark levels, as well as abso- 
lute sensitivity. Each of these corrections is performed by the 
IRAF^^ routine calacs (in stsdas.hst_calib.acs). 
So as to avoid the uncertainty that would otherwise have been 
introduced by the use of pre-launch flats in the OTFR pipeline, 
we updated the pfltfile header keyword for each raw image 
to specify that the newest available flat be utilized for each 
calacs run. 

IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy (AURA) under cooperative agreement with the National Science Foun- 
dation. 
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TABLE 1 
Observation Log 



Source 


2 


Obs. Date (UT) 


a (J2000.0) 


5 (J2000.0) 


Filters 


Line/Cont. 


A (A) 


Exp. Time (s) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


3C 33.0 


0.0597 


20 Aug 2006 


01 08 52.8 


+13 20 14 


FR647M 


Cont. 


5824 


1 x60 












FR716N 


Ha 


6950 


2x200 












FR551N 


[Olll] 


5302 


2x250 


3C 40.0 


0.0180 


17 Sep 2006 


01 20 34.0 


-01 20 34 


FR647M 


Cont. 


5593 


1 x60 










FR656N 


Ha 


6674 


2x200 












FR505N 


[Oiii] 


5092 


2x250 


3C 78.0 


0.0286 


13 Sep 2006 


03 08 26.2 


+04 06 39 


FR647M 


Cont. 


5654 


1 x60 










FR656N 


Ha 


6747 


2x200 












FR505N 


[Oiii] 


5147 


2x250 


3C 93.1 


0.2430 


8 Sep 2006 


03 48 46.9 


+33 53 15 


FR647M 


Cont. 


6842 


1 x60 












FR853N 


Ha 


8146 


2x200 












FR601N 


[Oiii] 


6228 


2x250 


3C 129.0 


0.0208 


25 Oct 2006 


04 49 09.1 


+45 00 39 


FR647M 


Cont. 


5615 


1 x60 












FR656N 


Ha 


6701 


2x200 












FR505N 


[Oiii] 


5112 


2x250 


3C 132.0 


0.2140 


20 Aug 2006 


04 56 43.0 


+22 49 22 


FR647M 


Cont. 


6677 


1 x60 












FR782N 


Ha 


7967 


2x200 












FR601N 


[Oiii] 


6078 


2x250 


3C 136.1 


0.0640 


4 Nov 2006 


05 16 03.1 


+24 58 25 


FR647M 


Cont. 


5852 


1 x60 












FR716N 


Ha 


6983 


2x200 












FR551N 


[Oiii] 


5327 


2x250 


3C 180.0 


0.2200 


11 Nov 2006 


07 27 04.5 


-02 04 42 


FR647M 


Cont. 


6710 


1 x60 












FR782N 


Ha 


8006 


2x200 












FR601N 


[Oiii] 


6108 


2x250 


3C 196.1 


0.1980 


3 Dec 2006 


08 15 27.8 


-03 08 27 


FR647M 


Cont. 


6588 


1 x60 












FR782N 


Ha 


7862 


2x200 












FR601N 


[OIII] 


5998 


2x250 


3C 197.1 


0.1280 


16 Sep 2006 


08 21 33.6 


+47 02 37 


FR647M 


Cont. 


6215 


1 x60 










FR716N 


Ha 


7416 


2x200 












FR551N 


[Oiii] 


5657 


2x250 


3C 219.0 


0.1744 


5 Dec 2006 


09 21 08.6 


+45 38 57 


FR647M 


Cont. 


6457 


1 x60 












FR782N 


Ha 


7704 


2x200 












FR601N 


[Oiii] 


5878 


2x250 


3C 227.0 


0.0858 


22 Jan 2007 


09 47 45.1 


+07 25 21 


FR647M 


Cont. 


5973 


1 x60 












FR716N 


Ha 


7127 


2x200 












FR551N 


[Oiii] 


5437 


2x250 


3C 234.0 


0.1848 


15 Dec 2006 


10 01 49.5 


+28 47 09 


FR647M 


Cont. 


6512 


1 x60 












FR782N 


Ha 


7770 


2x200 












FR601N 


[OIII] 


5928 


2x250 


3C 270.0 


0.0075 


25 Dec 2006 


12 19 23.2 


+05 49 31 


FR647M 


Cont. 


5541 


1 x60 












FR656N 


Ha 


6612 


2x200 












FR505N 


[Oiii] 


5044 


2x250 


3C 285.0 


0.0794 


11 Jan 2007 


13 21 17.8 


+42 35 15 


FR647M 


Cont. 


5934 


1 x60 












FR716N 


Ha 


7081 


2x200 












FR551N 


[Oiii] 


5402 


2x250 


3C 314.1 


0.1197 


30 Dec 2006 


15 10 22.5 


+70 45 52 


FR647M 


Cont. 


6154 


1 x60 












FR716N 


Ha 


7343 


2x200 












FR551N 


[Oiii] 


5602 


2x250 


3C 319.0 


0.1920 


12 Nov 2006 


15 24 05.5 


+54 28 15 


FR647M 


Cont. 


6556 


1 x60 












FR782N 


Ha 


7823 


2x200 












FR601N 


[Oiii] 


5968 


2x250 


3C 388.0 


0.0917 


29 Dec 2006 


18 44 02.4 


+45 33 30 


FR647M 


Cont. 


6004 


1 x60 












FR716N 


Ha 


7164 


2x200 












FR551N 


[Oiii] 


5466 


2x250 


3C 390.3 


0.0561 


30 Dec 2006 


18 42 09.0 


+79 46 17 


FR647M 


Cont. 


5808 


1 x60 












FR716N 


Ha 


6930 


2x200 












FR551N 


[Oiii] 


5287 


2x250 



Note. — The 19 3CR radio galaxies observed as part of the HST Cycle 15 SNAP program 10882 (PI: Sparks), listed by source name in ascending order (1) 3CR source name; (2) 
redshift; (3) observation date; (4) right ascension (J2000.0, in hours, minutes, and seconds); (5) declination (J2000.0, in degrees, arcminutes, and arcseconds); (6) ACS ramp filters used, 
listed for the 60 s continuum, the 400 s Ha, and the 500 s [Olll]A5007 observations, respectively; (7) Emission line (or continuum) observed, corresponding to filter configuration; (8) 
Redshift-adjusted Wavelength of emission line being observed, corresponding to the wavelength to which the ramp filter has been configured, with 2% bandpass for the Ha and [Olll] 
exposures, and 9% bandpass for the continuum exposures; (9) Exposure time in seconds. The continuum images are single exposures, while the emission line observations consist of 
two exposures in a two-point dither pattern. See §2 for more details on the sample selection and observations for this program. Also part of this program was 3C 371, the observation of 
which failed due to failure of guidestar aquisition. It is therefore excluded from presentation and analysis in this paper See Table[2]for a summary of the optical properties of the above 
sample, and Table[3]for a summary of its radio properties. 
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Fig. 3. — An illustrative example of the degree to which [Nil] flux varies relative to Ha among the SCR radio galaxies discussed in this paper. Here we 
plot de-reddened flux ratios of the diagnostic line [Nii]A6584 to Ha vs. redshift. The flux ratios are from the ground-based optical spectroscopy performed by 
Euttiglione et al. (2009) using the Telescopio Nazionale Galileo (TNG). The results from that work show that the relative contributions of [Nil] and Ha vary 
widely from source-to-source, and that [Nil] appears to dominate in most high-excitation galaxies (HEG, see Table |2). Red squares represent the four FR Is in 
our sample, while blue circles are FR lis and green triangles represent those objects exhibiting a CSS or unclassified radio morphology. 3C 132 is excluded from 
this plot as no TNG spectroscopy was available for this object. 
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Fig. 4. — A representative example of our iterative masking process for cleaning each image of cosmic rays and hot pixels. (Top left) Calibrated (by calacs, 
see §3.1) single-pointing Hq;-i-[Nii] exposure of 3C 33, one of two exposures that are combined during the mult idri z z le process to produce a single image (as 
described in §3.3). (Top right) One of the two CR and hot pixel masks created for the exposure presented at (Top left) by the cosmicrays routine, as discussed 
in §3.2. This mask was used to specify locations on the exposure in (Top left) at which to linearly interpolate across bad pixel regions by neighboring good pixels, 
as described in §3.2. In this representation, the bad pixels are displayed as black. (Bottom left) The same exposure as in (Top left) after two iterations of the 
CR/hot pixel cleaning process with the bad pixel masks. (Bottom right) The combined drizzled Hq;-i-[Nii] image of 3C 33, using the two dithered exposures (one 
of which is presented at Bottom left) that have been cleaned prior to the use of mult idri z z le by our iterative masking process. These images are significantly 
cleaner than those provided by the OTFR pipeline, particularly in the case of the single-pointing continuum images. See §3 for more information. 



In selecting the appropriate new flat to use (in lieu of the 
old ground flats used by the OTFR pipeline) we employ the 
same strategy that was used when the original ground flats 
for the ramp filters were created. That is, LRF flats are inter- 
polated from flats for ACS narrow, medium, and broad-band 
(i.e. non LRF) filters that correspond most closely with the 
redshift-adjusted wavelength of the emission line (or contin- 
uum) being imaged. The minor difference in central wave- 
length, throughput, and bandwidth of these (non LRF) filters 
corresponds to negligible differences in their corresponding 
flat fields (less than the ~ 7% uncertainty introduced in using 
the old ground LRF flats). 

With the replacement flats chosen on a case-by-case ba- 
sis and the image headers updated to reflect these changes, 
the calacs routine was run for each image, producing cal- 
ibrated FITS files (one for the continuum exposure, and two 
each for the Ha+[Nll] and [Olll] exposures) that have been 
flat-fielded and corrected for bias, dark current, and absolute 
sensitivity. It is upon these calibrated images that cosmic ray 
and hot pixel rejection was performed, as we discuss in the 
following section. 

3.2. Cosmic ray and hot pixel rejection 

The NOAO IRAF task cosmicrays 
(noao . imred . crut il) was run on copies of the science 
frames of each calibrated exposure (one for the continuum 
image, two each for the Ha+[Nll] and [Olll] images) to 
identify CRs and hot pixels. The routine identifies cosmic ray 
events by a detection algorithm based on a mean flux compar- 
ison within a specified detection window. A pixel within this 



window whose value is higher than a user-specified threshold 
dependent upon the mean value of its neighbors is considered 
"bad", and flagged for cleaning by the routine. Once the 
list of cosmic ray and hot pixel candidates is established 
over the whole image, cosmicrays replaces flagged pixels 
with the mean value of its four neighbors. A "cleaned" 
output file is then generated. We were careful to ensure that 
cosmicrays did not mis-identify real features (particularly 
galaxy nuclei or other physical, bright features) as CRs or hot 
pixels. The routine's parameters were adjusted if this was the 
case. 

The "cleaned" output from cosmicrays is insufficient 
for our needs, as the routine will often simply "dull" the CR 
rather than truly replace it. This is because cosmicrays 
flags pixels within a detection window that often contains 
neighboring pixels that are also affected by the CR (and 
are therefore far brighter than the true background). In the 
worst instances, depending on the parameters that are chosen, 
cosmicrays can sometimes reconstruct \hQ very CR it was 
attempting to mask. 

Therefore, we utilized the "cleaned" output from 
cosmicrays only as an intermediate step towards the 
creation of a bad pixel mask for each exposure. To that end, 
we divided the original, uncleaned copy of the calibrated 
science frame with the "cleaned" output from cosmicrays. 
The divided image was then examined by eye, and blinked 
with the uncleaned frame so as to qualitatively assess the suc- 
cess of cosmicrays in identifying offending pixels. The 
bad pixel mask was then created from the divided image, in 
which pixels whose values were above a specified threshold 
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Fig. 5. — A "before and after" example result from our continuum subtraction process, {a) Reduced FR716N //q;+[Nii] image of 3C 33, contaiminated by 
galaxy continuum, (b) The same image, with galaxy continuum subtracted via the method discussed in §3.5 



dictated those that were to be masked in the uncleaned image. 
This threshold was decided quaUtatively for each image 
based on our assessment of the success of cosmicrays 
in distinguishing truly bad pixels from real features. Often, 
accepting every flagged pixel from the divided image into 
the final mask would result in masking healthy pixels on the 
galaxy itself. If this appeared to be the case, the "allowance 
threshold" was set above the values of the incorrectly flagged 
pixels in the divided image, so that they were not made part 
of the final mask. We paid particular attention to flagged 
pixels within the immediate vicinity of the target galaxy, and 
ensured on a case-by-case basis that the pixels being masked 
were indeed CR events or hot pixels and not real features 
in the galaxy. This was accomplished through blinking 
individual exposures or by comparing the continuum image 
to archival if^Tdata from e.g. WFPC2 in the few cases where 
the correct choice was not obvious. 

The mask was then used with the f ixpix task to linearly 
interpolate values for masked pixels based on those of their 
nearest unmasked neighbors. While this is in effect similar 
to the procedure employed by cosmicrays, it is more ro- 
bust in that f ixpix is better able to replace a CR event with 
pixels generally consistent with that of the true background, 
cosmicrays, on the other hand, will often simply "dull" 
the CR rather than truly replace it, as it flags pixels within a 
detection window that often contains neighboring pixels af- 
fected by the CR (and are therefore far brighter than the true 
background). Moreover, the masking process allows the user 
to check at each step that no real feature on the galaxy is er- 
roneously altered. This process was repeated several times 
(typically three times) for each calibrated science frame (one 
continuum exposure, two Ha+[Nll] frames, and two [Olll] 
frames). At each iteration in the process we varied the detec- 
tion thresholds in cosmicrays, as well as the parameters 
associated with creating the subsequent mask, in order to en- 
sure that all CRs and hot pixels were ultimately identified and 
cleaned. In Fig.|4]we provide a representative example of the 
steps in our cosmic ray and hot pixel rejection strategy. The 
upper-left and lower-right panels in the figure serve as a "be- 
fore and after" comparison. 

In two cases (3C 33 and 3C 78) a large CR event landed 
directly on or near to the target galaxy image in the single 



pointing continuum exposure. In these cases we masked the 
CR while taking care to affect the galaxy light distribution as 
little as possible. Inevitably, some error due to pixel value 
interpolation remains on the final reduced images (as the con- 
tinuum exposure is used for both Ha+[Nll] and [Olll] in the 
continuum subtraction process, described in §3.4). However, 
while this error is difficult to characterize it is unlikely that 
it has a significant effect on the ultimate science value of the 
data, both for morphological studies of the emission-line re- 
gion and for e.g. flux measurements. 

3.3. Multidrizzle 

The calibrated, cleaned data were passed through the 
multidrizzle routine with the default parameters. The 
task calculates and subtracts a background sky value for 
each exposure, searches for additional bad pixels not already 
flagged in the data quality array, and drizzles the two ex- 
posures into outputs that are shifted and registered with re- 
spect to one another. From these drizzled exposures a me- 
dian image is created, which is then compared with original 
input images so as to reject the (very few) CRs not already 
cleaned in our previous CR rejection process, as we describe 
in the section above. More infor mation on the specifics o f 
multidrizzle can be found in iKoekemoer et all (l2002h . 
Final output images were left unrotated with respect to north 
to avoid pixel interpolation errors that might add uncertainty 
to our continuum subtraction and flux estimation strategies, 
which we describe in the sections below. 

3.4. Image Alignment 

To facilitate continuum subtraction it was first necessary to 
align the continuum exposures with their corresponding driz- 
zled emission line images to ensure registration down to a 
fraction of a pixel. A first pass registration was made by 
the world coordinate system (WCS) to achieve rough align- 
ment of the continuum image with its companion line emis- 
sion data. Subsequently, fractional shifts and rotations were 
manually applied using foreground stars and features on the 
galaxy as alignment aides (with the added help of the task 
imalign). Initial image registration was already well estab- 
lished for 8 of the 19 galaxies in our sample as their observa- 
tions were carried out using the same guide stars in all (con- 
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TABLE 2 
Optical Properties 





Ha Line Flux 


[Oiii] Line Flux 


log Lua 


log L[oiii] 




L.A.Shc 


Source 


(x 10 erg s ^ cm ^) 


(x 10 erg s ^ cm ^) 


f -1 \ 
(ergs 


(ergs 


Ionization Class ^' 


(kpc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


jj.yj 


0.428 


1.58 


Jy .J Z^J J 


40.0927 


HEG 


5.142 


3C 40.0 


1.69 


1.66 


39.0232 


39.0154 


LEG 


1.259 


3C 78.0 


2.18 


2.15 


39.5719 


39.5659 




2.208 


3C 93.1 


0.103 


0.070 


40.2147 


40.0470 




3.806 


3C 129.0 


0.12 


0.087 


38 0598 


37.9201 




0.243 


3C 132.0 


0.079 


0.081 


39.9738 


39.9847 


LEG 


9.77 


1 ^6 1 


0.205 


0.123 


39.2786 


Jy AjJKj I 




2.64 


1 Ro 


0.228 


0.488 


40 4605 


40.7910 


HEG 


16.78 


3C 196.1 


0.149 


0.136 


40.1756 


40.1360 


LEG 


8.72 


3C 197.1 


0.121 


0.091 


39.6925 


39.5687 


HEG 


3.04 


3C 219.0 


0.148 


0.143 


40.0509 


40.0360 


BLRG 


1.90 


3C 227.0 


1.931 


0.569 


40.5184 


39.9878 


HEG'^ 


21.76 


3C 234.0 


0.941 


2.050 


40.9068 


41.2449 


BLRG 


6.47 


3C 270.0 


1.91 


2.10 


38.3014 


38.3425 




0.63 


3C 285.0 


0.433 


0.711 


39.7925 


40.0079 


HEG 


6.49 


3C 314.1 


0.016 


0.070 


38.7325 


39.6781 




3.13 


3C 319.0 


0.027 


0.016 


39.4047 


39.4820 


LEG 


6.79 


3C 388.0 


0.056 


0.366 


39.0220 


40.1419 


LEG 


1.98 


3C 390.3 


11.5 


1.94 


40.9081 


40.4398 


BLRG 


2.97 



References. — fa^ lJackson & RawlingsHl997h . (b^ lButtidione et alJf2QQ^ . (c^ lPrietol fl993h 

Note. — (1) 3CR source name; (2) measured Ha (6563 A) line flux in erg s~^ cm~^; (3) total measured [Olll] (5007 A) line flux in erg s~^ cm~^; (4) total Ha luminosity in 
erg s~^ and (5) total [Olll] luminosity in erg s~^ (fluxes converted to luminosity using Hq -11 km s~^ Mpc~\ Q^m - 0.27, and Q.fy - 0.73); (6) Source classification of host (PR I = 
Fanaroff-Riley class 1 radio galaxy, QSO = quasar, BLRG = broad line radio galaxy, HEG = high-excitation FR II galaxy, LEG = low-excitation PR II galaxy, per the conventions in 
Qackson & Rawlings 1997); (7) largest measured angular size of line-emitting region (noted emission line in parentheses note the line for which the observed NLR appears largest, both 
lines win be listed if the NLRs cover approximately the same physical extent in both lines). 



tinuum, Ha+[Nll] , and [Olll]) exposures. The remaining 11 
observations included small slews of HST to allow for proper 
placement of the galaxy image on the WFC chip (necessary 
to image at the desired LRF wavelength). In these cases guide 
star acquisition was required between exposures. Registration 
for these 1 1 targets was therefore not initially ideal, and great 
care was taken to ensure proper alignment through manually 
applied shifts and rotations. After the appropriate shifts and 
rotations were applied, the WCS systems on the images were 
registered with respect to one another so as to enable proper 
WCS alignment when viewing the images. 

3.5. Continuum subtraction 

Here we describe the steps taken to produce Ha+[Nii] and 
[Olll] line emission images that are effectively free from con- 
tamination by continuum light from the galaxy. The aligned 
images were multiplied by their exposure times to set pixel 
values (in electrons per second for drizzled ACS images) to 
total electrons. The continuum image was then scaled by a 
factor determined for each image on a case-by-case basis. Ini- 
tially, this factor was roughly estimated using the synphot 
synthetic photometry package (specifically calcphot) in 
IRAK The calcphot models served as a qualitative guide 
in roughly establishing the scaling factor, which was later re- 
fined iteratively through a trial and error process that involved 
scaling and subtracting the continuum image from the line 
emission image, and examining residual pixel values on the 
resultant image in regions known to be dominated by con- 
tinuum emission on the original unsubtracted data. The pro- 
cess was repeated until these regions possessed a mode pixel 
value very close to the background level (which was effec- 
tively zero, after the multidrizzle process). Ultimately, the 
above strategies were employed in both qualitative and quan- 
titative ways on a case-by-case basis for each image, ensuring 
that the final subtracted image could be confidently consid- 



ered as "pure" line emission within the error established by 
the absolute photometry for this dataset. In Fig. [5] we provide 
a representaive "before and after" comparison in the case of 
continuum subtraction for 3C 33. 

3.6. Emission line luminosity measurements 

For each Ha+[Nll] and [Olll] image we have measured the 
total flux from all detected line emission in the galaxy's cen- 
tral regions. Emission-line fluxes were measured from the 
continuum- subtracted images using the apphot routines in 
IRAK We adopted an "all the flux you see" approach when 
establishing the size of the aperture with which the photomet- 
ric measurement was performed. In this way, we measure 
"total" flux as opposed to flux through an aperture that is the 
same across the the entire dataset. 

For consistency, we defined a "source region" by heavily 
smoothing a copy of each emission line image with a Gaus- 
sian kernel set to an arbitrarily large sigma. The aperture 
for the measurement was then set to the radius correspond- 
ing to a reasonable contour in the smoothed image enclos- 
ing all detected line emission from the galaxy. We also de- 
fined an outer boundary beyond which the image may be non- 
monochromatic (due to the ramp filters, see §2 for an expla- 
nation). The region between the "source region" and the outer 
boundary defined the area in which the Ha+[Nll] and [Olll] 
background level was measured. 

With these regions defined using the smoothed image, the 
corresponding regions on the original continuum- subtracted 
images were identified. The image was then examined to es- 
tablish the standard deviation of background, and the pixel 
values were thresholded to a value of order one sigma above 
the background level, corresponding to ~ 0.01 electrons per 
second. Counts above this threshold were considered signal 
and were counted in the photometric summation of flux val- 
ues, while pixels below this threshold were read as zero and 
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Fig. 6. — HST/ ACS narrow-band images of the Hq;-i-[Nii] (left panel, in red) and [Oiii] (right panel, in blue) line emitting regions in our sample of radio 
galaxies from the SCR catalog. Contribution from continuum has been subtracted using a strategy discussed in §3. In the bottom left corner of the left (right) 
panel we indicate the size of the field of view in arcsec (kpc). Both images are aligned, with East left and North up. The arrows in the upper right corner of some 
panels indicates the projected orientation of the radio jet axis on the sky, estimated from high resolution VLA radio maps (where available). Galaxies are listed 
in ascending order by SCR name. See Table[T]for more details on these observations, and §4 for descriptions of the emission line features of each target. 



were not counted. 

We summed counts above the pre-defined threshold within 
the source region (aperture), and then subtracted the esti- 
mated background measured similarly within the intermedi- 
ate region (exterior to the source region but interior to the 
outer boundary). Photometric conversion was then applied 
by scaling the residual (source minus background) value by 
the inverse sensitivity (the PHOTFLAM keyword), converting 
the value from electrons per second to flux units in erg cm~^ 
s~^ A~^ Flux measurements were made on unrotated images 
to avoid pixel interpolation errors. Emission line flux values 
were converted to luminosities via the luminosity distance. 
Non-rotated images were used for flux calibration to mini- 
mize pixel interpolation errors. Ultimately, the calibrations 
allow us to measure emission line flux with an uncertainty of 
- 15%. 

4. RESULTS & DISCUSSION 

In Figs. [6] -[8] we present the continuum-free line emission 
images listed by ascending 3C number. Ha+[Nll] images 
are shown to the left and are presented in orange and [Olll] 
images are shown to the right in blue. Where possible we 
have estimated the projected orientation of the radio jet axis 



on the sky from high spatial resolution archival Very Large 
Array (VLA) radio maps retrieved from the NASA/IPAC Ex- 
tragalactic Database (NED^^), indicated by an arrow in the 
top right corner of the galaxy's [Olll] image. This orientation 
was qualitatively estimated from the apparent inner jet axis 
for FR I radio galaxies, and along the hotspot-to-hotpsot axis 
for FR II radio galaxies. See Table [3] for a summary of the 
radio properties for objects in the sample. In instances where 
the figure does not display an estimate for the radio jet orien- 
tation, data was either not readily available, was of insufficient 
resolution to make a confident estimate, or was cospatial with 
the source on galactic scales (i.e. morphologically similar to 
CSS sources). 

4.1. Description of individual sources 

In the subsections below we discuss the observed morphol- 
ogy of the high surface brightness Ha+[Nll] and [Olll] line 
emission detected in each galaxy. Objects are listed in as- 
cending order by the name of the source. For information 
relating to the observations of each of the below galaxies, see 
Table [T] For quantitative data regarding the optical and radio 

http : / /nedwww . ipac . caltech . edu/ 
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properties of each object, see Tables [2] and [3l respectively. All 
images discussed below are presented in Figs. [6l [71 and[8l 

4.1.1. 3C 33; z = 0.059 

Extended regions of high surface brightness line emission 
are seen in both the Ha+[Nll] and [Olll] images. The low- 
and high-excitation regions trace a general "integral sign" 
shape extending ~ 5 kpc and oriented NE to SW. In both 
images a fainter, detached shell of low- and high-excitation 
emission is observed NE of the main "integral sign" feature. 
The [Olll] features more obviously clumpy, brighter hotspots 
along the lane of emission and near the nucleus than does 
the Ha+[Nll] image. Lon g slit spectrosc opy of the extended 
line emission obtained by lsimkini (Il979|) and lHeckman et al.l 
{1985) suggested that the line emitting gas is rotating along an 
axis oriented at a position angle (P. A.) of ~ 19° with respect 
to the FR II radio jet axis. (Fig.©. 

4.1.2. 3C 40; z = 0.011 

A dusty disk on 100 pc scales is seen in both Ha+[Nll] and 
[Olll] images, marked by the absence of emission along a 
circumnuclear rim that is more prominent in the [Olll] im- 
age than it is in Ha+[Nll] . Low- and high-excitation line 
emission is seen tracing the edges of the disk on the southern 
side, while [Olll] emission on the northern side of the disk 
appears less extensively distributed. The band marking the 
disk is not continuous, rather in both images it is interrupted 
by low-surface brightness line-emission apparently connect- 
ing the southern and northern sides of the disk in the same 
general location. Note that the apparent major axis of the 
disk appears roughly orthogonally oriented with respect to the 
FR II radio jet axis. (Fig.0). 

4.1.3. 3C78;z = 0.02S 

An optical jet is seen in both Ha+[Nll] and [Olll] emis- 
sion to extend from the core in a northeasterly direction (pro- 
jected on the sky). The bright, unresolved nucleus in both im- 
ages is surrounded by a generally isotropic region of narrow- 
line emitting gas in both Ha+[Nll] and [Olll]. iSparks et al.l 
feOOO) noted the presence of a face-on dusty nuclear disk in 
this FR I radio galaxy. Both distributions appear slightly more 
extended in the plane perp endicular to the opti c al synchrotron 
radio jet axis described by lSparks et al.l (Il995h . iPerlman et al.l 
(2006) notes the presence of an emission line knot cospatial 
with the radio jet axis. That we apparently detect both high- 
and low-excitation emission along the jet axis may be an im- 
portant result, and is worthy of future study. (Fig.©. 

4.1.4. 3C 93.1; z = 0.244 

A compact distribution of narrow-line emitting gas is seen 
surrounding the nucleus and extended ~ 3.8 kpc primarily to- 
wards the NE. The lopsided morphology observed in both the 
Ha+[Nll] and [Olll] emitting gas is evidence of the presence 
of an anisotropic nuclear ionizing radiation field. 3C 93.1 is a 
CSS source. (Fig.©. 

4.1.5. 3C129; z = 0.021 

Narrow-line emission associated with 3C 129 is only 
marginally detected in both filter configurations. Compact 
Ha+[Nll] emission with very low surface brightness is ob- 
served to be cospatial with the nucleus. Conversely, the nu- 
cleus is not detected in the [Olll] image, where instead only 
very faint extended emission that appears to be associated 
with the nucleus extends 1" to the north and south. 3C 129 is 
one of the four FR I radio galaxies in our sample (Fig. [6]). 



4.1.6. 3C132; z = 0.214 

A bright nucleus surrounded by fainter, more diffuse emis- 
sion is observed in both Ha+[Nll] and [Olll] images. The 
nucleus appears fainter in [Oiii] and what appears to be a 
faint "hotspot" is located directly to the north of the nucleus. 
(Fig.©. 

4.1.7. 3C136.1; z = 0.064 

In Ha+[Nll] , an elongated nucleus is seen with major axis 
oriented north- south on the sky. Faint, disjoint distributions 
of narrow-line gas surround the nucleus in two filamentary 
plumes extending from the core along a N-S axis roughly per- 
pendicular to the FR II radio jet axis. Clumpy patches of dis- 
parate, low-surface brightness emission are observed at radii 
of ~ 1''. In [Olll] there is an apparent double nucleus with 
hotspots more clearly separated than in the Ha+[Nll] image. 
The north- south filaments seen in HQf+[Nii] are not detected 
in [Olll]. (Fig. El). 

4.1.8. 3C180; z = 0.220 

A bi-conical series of shells is observed in both 
Ha+[Nll] and [Olll], reminiscent of other bi-conical NLR 
morphologies observed in some Seyfert II g alaxies (e.g., Mrk 
3, Mrk 573, NGC 5252, ICapetti et al.l [T999). At both low- 
and high-excitation we observe a bright, lopsided emission 
region cospatial with the nucleus. Fainter "plume-like" fea- 
tures extend from both the northern and southern extremities 
of the extended core of emission. The bi-conical series of 
shells of emission extends ~ 8 kpc to both the NE and SW 
of the core. The abundance of optical line emitting gas along 
a preferred axis suggests the presence of a highly anisotropic 
ionizing radiation field. 3C 180 lacks a Fanaroff-Riley clas- 
sification, and we are unable to characterize the morphology 
of its radio sour ce in a ny meaningful way given the available 
data. McCarth v et al.1 (fl995.) note that 3C 180 resides in a 
cluster. (Fig.|7]). 

4.1.9. 3C196.1; z = 0.m 

A high surface brightness lane of Ha+[Nll] emission 
"snakes" NE from the bright elongated core, reminiscent of a 
"tadpole" shape. Ha+[Nll] emission is also observed extend- 
ing from the core along an axis perpendicular to the axis of 
elongation of the "tadpole"-shaped distribution of ionized gas. 
In the [Olll], high excitation emission is observed in two lo- 
calized hotspots cospatial with the elongated bright core seen 
in Ha+[Nll] , and the "tadpole" structure is not detected, nor 
is the emission along the perpendicular axis. (Fig.|7]). 

4.1.10. 3C197.1; z = 0.13 

Bright, unresolved nuclei are seen in both images, sur- 
rounded on all sides by a fainter distribution of line-emitting 
gas. The nucleus appears larger in Ha+[Nll] than it does in 
[Olll]. (Fig.Cl). 

4.1.11. 3C 219; z = 0.174 

Compact emission is seen in both bands, as is faint, dif- 
fuse emission surrounding the unresolved source at the center. 
(Fig. El). 

4.1.12. 3C 227; z = 0.0^6 

Compact emission cospatial with the galaxy's nucleus is 
seen at both low- and high- excitation. In [Olll] we observe 
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an extremely large arm of emission extending ~ 20 kpc, rep- 
resenting by far the largest distribution of line-emitting gas 
in our sample. The large structure is oriented along an axis 
roughly perpendicular to that of the FR II radio jet. .Prieta 
(11993) studies this object in great detail. (Fig. [7]). 



nucleus from both sides of the disk, and are elongated along 
the direction of the jet (east- west on the sky, and nearly per- 
pendicular to the major axis of the dusty disk). Nuclear [Olll] 
emission is also seen extending from both sides of the disk, 
though seems to be largely absent on the disk itself. (Fig. [8]). 



4.1.13. 3C234; z = 0.lM 

The prominent tidal arm described bv lCarletonetaP (119841) 
is seen in high surface brightness Ha+[Nll] and [Olll] emis- 
sion. To the east is the bright quasar-like nucleus of the 
galaxy, and in both images a bright feature is seen extend- 
ing 0.5'' from its northwestern edge. The overall elongation 
of the emission line region, including the tidal arm, is in a 
projected orientation that is roughly parallel with that of the 
radio jet axis on large scales. (Fig. [8]). 

4.1.14. 3C 270; z = 0.0077 

The ~ 120 pc dusty disk originally studied by iJaffe et al.l 
(1 19931) is notably absent of emission in both Ha+[Nll] and 
[Olll], though there remains some Ha+[Nll] emission along 
the inner regions of the disk. The disk is largely edge-on 
with respect to the line of sight though is inclined such that 
the western side of the disk "faces" the observer slightly. 
Note that the western half of the galaxy is noticeably brighter 
than the eastern half. Cones of high surface brightness 
Ha+[Nll] emission are seen extending from the unresolved 



4.1.15. 3C 285; z = 0.079 

A complex, filamentary system of dust lanes bisects the nu- 
clear region of the galaxy in a northeasterly direction and per- 
pendicular to the line of sight. The extinction associated with 
the predominant lane in this system is seen in the Ha+[Nll] , 
where we observe an extended region of high surface bright- 
ness emission with complex morphology. A bright, dense col- 
umn of line-emission extends ~ 1'' northward from the nu- 
cleus at the center of the image. Slightly to the southeast a 
small, bright "integral sign" feature is seen in general align- 
ment with the large dense column. Both of these features are 
surrounded by fainter narrow-line emitting gas arranged in fil- 
amentary plumes as well as clumpy patches. In the [Olll] im- 
age we observe a dramatically different morphology than is 
apparent from the Ha+[Nll] . A high surface brightness lane 
of emission "snakes" westward from the center of the galaxy 
to form a disjoint "W" shape. The nucleus of the galaxy, such 
as it is, is not discernible from the image. Bright patches of 
emission populate the prominent filament at regular intervals 
along its length. A region of particularly high surface bright- 
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ness emission is noted at its westernmost extremity. (Fig. [8]). 
4.1.16. 3C 314.1; z = 0.ll9 



Diffuse Ha+[Nll] emission is distributed in a complex mor- 
phology, while the [Olll] emission is decidedly more com- 
pact, albeit with faint, wispy plumes extending from the nu- 
cleus and slightly toward the southwest. (Fig. [5]). 

4.1.17. 3C 319; z = 0.192 

Very faint line emission is detected in both images, dis- 
tributed in a diffuse, complex morphology at extended dis- 
tances from the nucleus. The [Olll] is decidedly more com- 
pact than is the Ha+[Nll] . (Fig. [5]). 



4.1.18. 3C 388; z = 0.092 

A bright unresolved nucleus is seen in both images. Small, 
low surface brightness "plumes" can be seen extending to- 
ward a northeasterly direction in the Ha+[Nll] , while a more 
prominent conical feature is seen on the southwestern side of 
the nucleus in the [Olll] image. (Fig. [5]). 

4.1.19. 3C 390.3; z = 0.056 

Very bright Ha+[Nll] emission dominates the image from 
the unresolved nucleus. The noticeably rhomboidal shape of 
the nucleus in the [Olll] is likely a result of the diffraction 
spikes seen in both images. (Fig. [8]). 

4.2. Discussion: correlations 
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vice-versa. For a discussion on the implications of these emission line luminosity ratios, see §4.2. 



Below we discuss significant correlations, consistent with 
the results of previous related studies, that are apparent in our 
data. In §4.3 we discuss the physical interpretation of these 
observed trends, and compare our results to past literature. 

4.2.1. Emission line luminosities vs. redshift 

In Fig.[9]2 and [3? (top left and top right, respectively) we 
plot the measured Ha+[Nll] and [Olll] luminosities {Lwcx and 
^[Oiiii) vs. redshift for the galaxies in our sample. FR I radio 
galaxies are plotted in red squares, while FR lis are plotted in 
blue circles and the three galaxies whose radio morphologies 
are uncategorized are plotted in green triangles (this conven- 
tion holds throughout this paper). As expected, the emission 
line luminosities follow a positive redshift-luminosity trend. 



albeit with a greater degree of scatter than that found in the 
canonical K-z relation for radio galaxies (where K is the 
NI R ^-band magnitu de of the elliptical radio galaxy host, 
see lWillott et al .112003 for a detailed study of this correlation). 
This higher degree of scatter is not necessarily surprising for 
a number of reasons, some of which we summarize below: 

1 . Warm optical line emitting gas in radio galaxies is nat- 
urally expected to inhabit a greater range in luminos- 
ity space than are the NIR starlight distributions of 
their elliptical hosts. As the stellar component of ra- 
dio galaxies is believed to have formed at high redshift 
and evolved passively in the time since (barring other 
effects), a tight clustering of ^-band magnitudes is ex- 
pected. The amount, ionization state, and detectability 
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Fig. 10. — A representative comparison of data quality from the WFPC2 LRF emission line survey of SCR radio galaxies as presented inlP rivon et al. ! (120081) 
in comparison with the data presented in this paper, (a) Continuum- subtracted WFPC2 FR533N 2 x 300 s exposure of the [Oiii] emission in 3C 234; (b) 
continuum-subtracted ACS WFC FR601N exposure of the [Oiii] emission in 3C 234, from the dataset presented in this paper. See Table|2]and §2 for details. 



of gas in the nuclei of radio galaxies, however, is sub- 
ject to a number of factors that greatly influence overall 
emission line luminosities. 

2. Our sample includes low-power FR I radio galaxies, 
characterized by narrow-line emission, and FR II radio 
galaxies and QSO characterized by broad optical line 
emission. Broad line radio galaxies (BLRGs) generally 
have an optical excess, resulting in a larger contribu- 
tion to emission line flux. Not surprisingly, the BLRGs 
in our sample (3C 33, 197.1, 219, 285, 234, & 390.3) 
generally inhabit the upper regions of Fig. ^ and [2?, 
contributing to the overall scatter. 

3. These galaxies range in redshift from 0.0075 < z < 
0.224 and were observed with a variety of LRF filter 
configurations (see Table [2]). This results in varying 
amounts of emission line flux loss beyond the edge of 
the passband, as well as contamination from [Nil] . 

In Fig. ^ we plot our Ha+[Nll] luminosities in blue 
squares vs. redshift, overplotted with the measured and ex- 
trapolated Ha+[Nll] luminosities for the sample of 3CR ra- 
dio galaxies from Privon et al.l (I2008D . The slightly larger 
scatter observed in the Privon data reflects the fact that 
roughly 50% of their Ha+[Nll] luminosities have been esti- 
mated from measured [Olll] luminosities using ratios found 
in the literature. In general, however, our galaxies follow a 
redshift-luminosity trend roughly consistent with that found 
for the (generally) higher redshift subset of the 3CR studied 
in Privon et al. (2008). See Fig. [TO]for a representative and 
qualitative comparison of WFPC2 and ACS data quality from 
the Privon observations and this paper, respectively. 

4.2.2. Emission line ratios 

In Fig. [9]i (bottom right) we plot the measured [Olll] 
luminosities vs. the measured Ha+[Nll] luminosities for 
the objects in our sample. Those galaxies lying above the 
black diagonal line (tracing L[oiii] = ^Ha) possess [Olll] lumi- 
nosities greater than their Ha+[Nll] luminosities (and vice- 
versa). Note that, in general, Lho and L[oiii] are roughly sim- 



ilar across the sample, with some notable exceptions whose 
names we have indicated in the figure. 

In future papers it will be necessary to investigate emission 
line ratio maps, which may carry important implications with 
regards to the orient ation-base d unification scheme of radio- 
loud AGN (e.g., Pad ovani & U rrv 1992; Chiaberge 2004, and 
references therein). It is expected that low-excitation lines 
(like Ha+[Nll] ) are emitted further from the AGN than are 
higher excitation lines like [Olll], as the fraction of highly ex- 
cited gas will depend on the proximity of that ga s to the ioniz- 
ing source, in this case the AGN (e.g., Jackson ^ Rawlingsl 
1997, and references therein). This spatial dependence of 
^[0111] /^Ha also leads to a dependence upon the sources of 
obscuration present in the galaxy. In the unification model 
for radio galaxies, wherein FR Is are unified with BL-Lacs 
and FR lis are unified as the parent population of steep spec- 
trum QSO, we expect low-excitation gas to lie further from 
the obscuring torus, while the fraction of high-excitation gas 
observed will follow a tighter dependence upon the angle at 
which the observer views the obscuring torus. We may there- 
fore expect that intrinsically similar radio-loud AGN viewed 
at a variety of angles will have unequal [Olll] luminosities 
but roughly equal H a+[Nll] luminosities (Baum & Hec kmanI 
ll989bl:IJackson & Br owne 1990). 

It will also be important to investigate whether [Olll]/Ha 
ratios of order unity are indicative of Seyfert-like nebular ac- 
tivity. We note that, in general, the majority of objects in our 
sample with L[oiii] /^Ha ^ 1 in (i.e. those galaxies significantly 
offset from the diagonal line in Fig.|9]i) are BLRGs, with the 
exception of one FR I. Though we are unable to make statis- 
tically significant inferences with respect to these issues (as 
our sample consists of only 19 galaxies), further investigation 
into these line ratios and associated line ratio maps may yield 
important results. 

4.2.3. Emission line luminosities vs. total radio power 

In Figs. [TlT a) and fTlT /?) we plot emission line luminosity 
vs. total radio power at 178 MHz for Ha+[Nll] and [Olll], 
respectively. We note a clear upward trend in both figures, 
consistent w i th the results of the comprehensive study by 
IWillott et"an (Il999h . which also found a tight correlation be- 
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tween emission line and radio luminosity in flux limited sam- 
ples of radio galaxies. 

We discuss interpretations of this result in §4.3, but it is first 
important to frame the results of the above emission line lumi- 
nosity comparisons in the context of ionized gas distribution 
morphologies. 

4.2.4. Morphology of the extended narrow line regions 

The low- and high- excitation narrow-line regions we ob- 
serve in our sample, overall, exhibit extended and complex 
morphologies characterized by (e.g.) clumpy regions of emis- 
sion (e.g. 3C 33), whispy tendrils (3C 136.1, 3C 319), ex- 
tended filaments and lanes (3C 196.1, 3C 234), and bi-conical 
series of shells (3C 180). In general the morphologies ob- 
served in [Olll] mirror those in Ha+[Nll] with some impor- 
tant exceptions (e.g. 3C 136.1, 3C 196.1, 3C 227, 3C 285). 
We also note the presence of bright unresolved nuclear emis- 
sion in the cores of every object with the exception of 3C 129 
in [Olll] and 3C 285 in Ha+[Nll] and [Olll] (no clearly "point 
source-like" nucleus is apparent in our images for that object). 

In Figs. [T2ta) and[T2tZ?) we compare the projected linear 
sizes of the line-emitting regions with both redshift and el- 
liptical host galaxy effective radius (respectively). We have 
estimated the gas distribution sizes, in kpc, from the largest 
angular extent (L.A.S.) subtended by the high surface bright- 
ness emission (as projected on the sky). While this is a rough 
estimate, it provides important insights as evidenced by the 
weak upward trend of L.A.S. with redshift (Fig.[T2]2), imply- 
ing that larger distributions of optical line emitting gas are 
found at higher redshifts. 

Values for the e ffect ive radii 7?eff in Fi^. \Vl\ b) are from 
iFlovd et al.l (l20Q8h and iDonzeUi et aP (l2007h . and were de- 
rived from surface brightness profile fitting to if^T/NICMOS 
//-band images (and some NIR imaging from TNG, in the 
case of Donz elli et al.ll2QQ7i) . While the methods employed 
by these two studies differ, their results are generally consis- 
tent, and any uncertainty introduced in plotting results from 
both papers is reflected in the characteristic error bars we have 
indicated in the figure. In light of these added uncertainties, 
however, we are unable to make any significant inferences, 
though L.A.S. and R^a do appear weakly correlated for galax- 
ies with 7?eff < 4 kpc, while no such correlation is evident for 
those objects with 7?eff > 4. 

4.3. Physical interpretation of the trends observed 
We have discussed three trends apparent in our data: 

1. Emission line luminosity (both for Ha+[Nll] and 
[Olll]) appears positively correlated with the redshift 
of the host galaxy. 

2. The Ha+[Nll] and [Olll] luminosities are also tightly 
correlated with total radio power at 178 MHz. 

3. Projected sizes of the ionized gas distributions in our 
sample are positively correlated with host galaxy red- 
shifts. That is, we observe larger line-emitting regions 
at greater distances. 

The abov e resul t s are consistent w it h tho se of (e.g.) 
IWiUott et al.l ([T999I [2003h : iPrivon et al.l (l2Q()8h (and refer- 
ences therein), as well as many other past studies of emis- 
sion line gas in radio galaxies. The tight redshift-luminosity 
correlation in flux-limited samples results in the selection 
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33.65 


3C 40.0 
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32.20 


3C 78.0 
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15.0 
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3C 93.1 
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14.5 
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9.5 


33.55 


3C 219.0 


FR II 


44.0 


34.54 


3C 227.0 


FR II 


28.0 


33.74 


3C 234.0 


FR II 


29.0 


34.47 


3C 270.0 


FRI 


44.0 


32.75 


3C 285.0 


FR II 


10.5 


33.23 


3C 314.1 


FRI 


9.0 


33.59 


3C 319.0 


FR II 


17.0 


34.20 


3C 388.0 


FR II 


22.0 


33.69 


3C 390.3 


FR II 


44.0 


33.55 



References. — Fanaroff & Rilev 1974; Spinrad et al. 1985 

Note. — (1) Source name; (2) Fanaroff-Riley classification; (3) 178 MHz flux 
density in Jy; (4) Logio of radio power at 178 MHz, in erg s~^ 

of intrinsically brighter objects at higher redshifts. In ra- 
dio flux density selected samples such as the SCR, this is 
manifest in more powerful radio galaxies (FR lis) inhabit- 
ing a higher redshift range than do their intrinsically lower- 
power counterparts (FR Is). Moreover, the alignment ef- 
fect discussed in §1 suggests that propagation of the radio 
jet pl ays a role in the excitation of the ISM into line emis- 
sion (iBaum & HeckmanI Il989al: iMcCarthvl 1 19931: iBest et al.1 
It is therefore unsurprising that more powerful ra- 
dio galaxies appear associated with higher emission line lu- 
minosities, as they are able to collisionally excite and pho- 
toionize greater amounts of ga s than a re lower -power (an d 
therefore lower-redshift) AGN (iBaum & Heckm^ fT989allbl : 
iRawlings & Saundersl ll99ll:llnskip eTa l. 2002a,lJ7^As accre- 
tion disk radiative l uminosity appears r elated to t he jet kinetic 
energy flux (iBaum & Heckmanlll989bl: Rawlings & Saunderi 
19911 this can be used to place constraints on the physics of 
the AGN. The greater mechanical energy input and number 
of ionizing photons provided by more powerful radio sources 
also explains the dependence of emission line region size on 
redshift: we expect the more powerful radio galaxies at higher 
redshifts to be capable of exciting gas at greater distances 
from the central engine, resulting in the observed trend. 

It should also be noted that dense clouds of gas are ob- 
served on increasingly larger scales at higher redshifts, ulti- 
mately forming the vast Ly-a hal os that surr ound the most 
powe rful radio galaxies ( van Qji k et al.l 1 19971: [Reuland et al.l 
l2QQ3h . While this dependence of gas cloud density and extent 
on redshift is important for objects with z > 0.6 (the point 
near which the alignment effect precipitously tightens), it is 
not likely to contribute in any meaningful way to our compar- 
atively low redshift sample at z < 0.3. 

5. SUMMARY AND CONCLUDING REMARKS 

In this paper we have presented HST/ACS narrow-band 
imaging of Ha+[Nll] and [Olll]A5007 line-emission in the 
central regions of 19 low-redshift radio galaxies from the 3CR 
catalog. We have chosen the 3CR as the basis for our sample 
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Fig. 1 1 . — Emission line luminosity vs. total radio power at 178 MHz for (a) Hq;+[Nii] and (b) [Oiii]. Radio luminosities are from'Spinrad et al." fl985). Red 
squares represent the four FR Is in our sample, while blue circles are FR lis and green triangles represent those objects exhibiting a CSS or unclassified radio 
morphology. We note the clear upward trend in both figures, consistent with the results of the comprehensive study by Willott et al. (1999) which found a tight 
correlation between emission line and radio luminosity in flux limited samples of radio galaxies. See Table[3]for a summary of the radio properties in our sample. 
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Fig. 12. — Comparison of the projected linear sizes of the extended high surface brightness Ha+[Nii] emission line regions with both id) redshift and ib) 
effective radii R^a of the host galaxies. The size estimates represent the largest angular size (L.A.S) of the narrow line regions detected in our Hq;+[Nii] imaging. 
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their R^ff was available. Note the general upward trend of the L.A.S. with redshift, as well as a weaker upward trend with /?eff- The former result is consistent with 
the emission line study of Privon et al. ( 2008), and is expected given the tight redshift-luminosity correlation resulting in the selection of higher power sources at 
higher redshifts. The higher-power sources are likely able to ionize gas at greater distances from the central engine, contributing to the upward trend we observe. 



as it is radio flux density selected at a frequency dominated by 
unbeamed radio lobes, and is therefore free from orientation 
bias with respect to all //^T wavelengths. Moreover, the SCR 
has been extensively covered by past ground- and space-based 
imaging and spectroscopy programs, yielding a robust, cross- 
spectrum database containing galaxies exhibiting a wide va- 
riety of intrinsic characteristics. We limit our sample to the 
low-redshift subset {z < 0.3) of the SCR to maximize spatial 
resolution, as well as to ensure that the emission lines fall 
within the redshift range of high sensitivity. The ACS obser- 



vations were carried out during HST Cycle 15, and in total, 
19 galaxies were observed prior to failure of the instrument in 
January 2007. Reduction of these data was performed inde- 
pendently of the automated OTFR pipeline so as to improve 
on flat-fielding and to more proactively clean the images of 
cosmic rays and hot pixels. In addition, contribution from 
continuum was subtracted to yield effectively pure emission 
line images. 

We have compared the resulting emission line luminosities 
with both redshift and total radio power, finding weak positive 
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correlations for each. We have also noted an upward trend be- 
tween the projected size of the emission line gas distributions 
and redshift. These results are consistent with those from past 
studies, and we have interpreted the apparent trends to be a 
result of higher power radio galaxies, preferentially found at 
greater distances, which can shock and photoionize greater 
amounts of the ISM further away from the nucleus, thus re- 
sulting in brighter and more extended emission-line regions. 

Regardless of the interpretation, recombination times 
for warm gas {T ^ 10"^ K) are of order 10^ years 
(lOsterbrock & F erlandl2006). while lifetimes of radio sources 
are of order 10^ years. The process by which the gas is ion- 
ized should therefore be ongoing throughout the lifetime of 
the source, implying a strong connection between AGN activ- 
ity and observed emission line properties. The sensitive, high 
spatial resolution images presented in this paper may there- 
fore enrich future works studying jets propagating through 
the ISM, and their relationship with extended and compact 



star formation regions as well as X-ray coronal halos. Future 
emission line studies of larger samples of radio galaxies may 
also deepen understanding of radio-loud unification schemes, 
as warm optical line emitting gas traces fundamental energy 
transport processes coupled to feedback at the interface of the 
AGN and its surrounding medium. 
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